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Abstract 

The Wigner functions on the one dimensional lattice are studied. Contrary 
to the previous claim in literature, Wigner functions exist on the lattice with 
any number of sites, whether it is even or odd. There are infinitely many 
solutions satisfying the conditions which reasonable Wigner functions should 
respect. After presenting a heuristic method to obtain Wigner functions, 
we give the general form of the solutions. Quantum mechanical expectation 
values in terms of Wigner functions are also discussed. 

PACS: 01.55.+b; 03.65.-w; 03.65.Sq 
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I. INTRODUCTION 



The Wigner function, the Wigner transform of the density matrix, has been recognized 
as a bridge between quantum and classical mechanics It is the quantum mechanical 
counterpart to the distribution function in the classical phase space: the Wigner function 
W(q,p) gives the position distribution if integrated over the momentum p and the momen- 
tum distribution if integrated over the coordinate q. Though the Wigner function may be 
negative, the equations of motion for W(q,p), in the classical limit, are reduced to the clas- 
sical Liouville equation for the distribution function. The Wigner functions have been used 
to analyze various phenomena where the quantum and classical correspondence is crucial 
|J. One of the important and interesting applications of the Wigner functions is the 
quantum tomography, where the Wigner function W(q,p) for a certain quantum system is 
measured for any q and p and thus the full quantum information of the system, density 
matrix, can be reproduced || p. 

Some time ago Wootters and Cohendet et. al. |7| studied Wigner functions for a finite 
discrete space like a spin space. For an integer spin s in a pure state, they found a possible 
Wigner function is given by 

1 s ( Att \ 

W ( a > & ) = 2s + 1 ^ 6XP V 2s + 1 / + ( a ' & = _s ' _s + X ' • • • ' s )' W 

where i/j(a) is the amplitude for each spin state labeled by an integer a = —s, — s + 1, . . . , s 
and ip(a + 2s + 1) = ip(a) is assumed. Their results, however, can not be applied to half- 
integer spin systems. The total number of sites is N = 2s + 1 and the Wigner function given 
in || and J7J is only for the lattice with odd number sites. 

Leonhardt proposed the Wigner function for an even number of sites N and discussed 
applications to the quantum tomography in a finite discrete space ||. However, the Wigner 
function given in || involves "ghost variables" with vanishing probability. Thus the con- 
struction of Wigner functions on a discrete space or lattice has still been an open problem. 
The relations between various quasidistribution functions including Wigner functions and 
the finite discrete space with arbitrary Floquet angles have also been discussed [[J fLOf . 



In this paper, we relabel the spin state a by an integer value a = s + a with a ranging 
from to 2s = N — 1 (the total number of sites is iV = 2s + 1). The equivalence between 
these two labeling will be discussed in Sec. 5 in detail. Labeling each site by an integer 
(0, 1, . . . , N—l), the Wigner function of Wootters and Cohendet et. al. is given by 

1 N ^ / 4?r \ 

W ( a > h ) = E ex P \ % ^ hh ) <^ a + h M a -h), (a,b = 0,l,---,N-l), (2) 

where 0(a + N) = 0(a) is assumed. Generally Wigner functions are written in terms of the 
Fano operators (matrices) A (a, 6) J l I 



as 



I N-l 

W(a,b) = - 0>i)A aia2 (a,&)0(a 2 ), (a, b = 0, 1, • ■ • , N - 1), (3) 

iV aia 2 =0 



or in terms of the density matrix p as 
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W(a,b) = ^TrpA(a,b). (4) 

In this paper we postulate the following four conditions (A)~(D) for eligible Wigner 
functions and corresponding Fano operators. 

(A) The Wigner functions reproduce the distribution in the configuration space: 

N-l 

£ W(a,6) = |0(a)| 2 , (5) 

6=0 



or equivalently 



N-l 

E A a ia2 (a, b) = N5 a>ai 5 ai>a2 . (6) 

6=0 



(B) The Wigner functions reproduce the distribution in the momentum space: 



N-l 



£ W(a,6) = |0(6)| 2 , (7) 



a=0 



where the Fourier transform <fi is given by 

<K & ) = E ex P (- i ^ ba ) <^ a )> (6 = 0, 1, . . . , TV - 1), (8) 



or equivalently 



.Z7T 



a=0 



J2 A aia2 (a,6) = exp [i— b(a x - a 2 ) ) , (9) 



(C) The Wigner functions are real. Therefore the Fano operators are hermitian: 

W(a, by = W(a, 6), or A(a, 6)+ = A(a, 6). (10) 

(D) The Fano operators are orthonormalized and complete: 

Tr(A + (a,6)A(a',6')) = N5 aa ,5 bb , (11) 

The completeness in the last condition (D) is necessary for that the density matrix can be 
expressed in terms of the Wigner functions. The completeness, together with the orthonor- 
mality, enables us to obtain the simple form: 

JV-l 

p ai ,a 2 = E W(a,b)A aia2 (a,b). (12) 

a,fe=0 



3 



We note that Eqs.([|) and (|T2| ) are valid for a mixed state as well. It is easy to see that 
these four conditions (A)~(D) are satisfied by the Wigner function given by Wootters and 
Cohendet et. al. 

In this paper we will show that the Wigner functions with the properties (A)~(D) exist 
for the lattice with any number of sites, and therefore any spin system. We do not need any 
"ghost variables" as in ||. In fact there are infinitely many solutions and we will give the 
general form. In Sec. 2 we will first give a heuristic way to obtain a solution for even site 
systems. The general form for an arbitrary number of sites will be given in Sec. 3. In Sec. 4 
we will discuss expectation values of physical quantities in terms of Wigner functions and 
the problem on the two ways of site labeling. 

II. CONSTRUCTION OF A WIGNER FUNCTION ON AN EVEN SITE LATTICE 

In this section, we describe a heuristic method to construct a Winger function on an even 
site lattice and show its existence. Let N be an even integer and Z N be the cyclic group 
whose elements are labeled from to N — 1. We construct the N x N Fano operators by 
a reduction method starting from 2N x 2N matrices. The enlarged matrices have a similar 
form to a cyclic representation of the Weyl matrices on an odd site lattice. 

We define phase and shift operators P and S as 

and two operators T and B as 

T ai a 2 = $a 1 ,2N-l-a 2 i B aia2 = <5 aiia2 + (_i)a 2 , (14) 

which are 2N x 2N matrices, where ai, a 2 £ Z 2 n and u' is the primitive 2A-th root of unity: 

^ = ex P ■ ( 15 ) 



V2A, 

These operators have the following properties: 

p + = p -i = p * ; s+ = s -i = s t^ T+ = T -i = T = T * ; 

PS = lj'- 1 SP, TS = S^T, BT = TB. (16) 

We consider a set of matrices W(a, b) constructed from these matrices: 

W(a,b) =u'- ab P b S- {2a+l \ 
W(a + A, b) = J-^+N)bpb s -{2a+i) B ^ ( 17 ) 

where a £ Zat,6 £ They are orthogonal in the trace norm and complete in the 

2A x 2A-dimensional operator space, but not hermitian. We define matrices A(a, b) by 

A(a, b) = W(a, b)T, (18) 

where a, b £ Z 2 jy. For these operators A (a, b), it can be shown that the condition (D) in the 
previous section is fulfilled and the following properties (A'), (B') hold. 
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(A') A property which is similar to (A) in the previous section: 



2JV-1 



A aia2 (a, b) = N5 aai 5 aia2 , 
A aia2 (a + N,b) = N5 aai 5 aia2 , (19) 



b=0 

2N-1 



6=0 

where a G Zjy. 



(B') A property which is similar to (B) in the previous section: 

2N-1_ / 27T \ 

^ A aia2 (a, b) = exp ( i— b(a x - a 2 ) ) , for even b 

a=0 ^ ^ V ' 

2N-1 

]T A aia2 (a,6) = 0, for odd b (20) 



a=0 

where b G Z 2N . 

Note that the Weyl matrices on an odd lattice, i.e. for odd N, can be represented as 

W(a, b) = LU -^p2b s -2a = s - ap 2b S -a^ ^ 

using similarly defined N x N phase and shift operators, where u is the primitive iV-th root 
of unity. 

We restrict the complex 2iV-dimensional vector space, on which the A(a, b) matrices act, 
to the iV-dimensional £> -invariant subspace V: 

V = {(j) G C 2N : B<P = 0}, (22) 

which means 0(2a) = 0(2a + l),a G Z N if G V. The matrices A(a, b) take N x N- 
dimensional form on this restricted space, e.g. choosing (j>(2a), a G Z N as independent ones, 
and we denote the reduced matrices by A'(a,b). The hermiticity is recovered after the 
reduction. It can be seen that the next two identities hold for the A', 

A'(a,b) = (-l) b A'{a + N,b), (23) 

(A'(2a, b), a G Z N , b G Z 2N ) = (A'(2a + 1, b), a G Z N , b G Z 2N ), (24) 

where ( ) means the vector space spanned by the elements between the bra and ket. Only 
N 2 matrices out of AN 2 matrices of A' are independent. We define Fano operators as 

A(a, b) = \ {A'(2a, 2b) + A'(2a, 2b + 1)} , (25) 

where a, b G Z N . We can see that these operators satisfy the properties (A)~(D). The 
explicit form of the corresponding Wigner functions is given by 

1 JV-l 

W(a, b) = ^-J2 ^' 4kb + u' 2h{2b+1) W(a + h)<p(a - h) 

1 JV-l 

+ ^7 E {u' {2h+1)2b + u'( 2h+1 ^ 2b+1 ^*(a + h)<P{a -h-l). (26) 
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III. GENERAL FORM OF WIGNER FUNCTIONS ON A LATTICE 



In this section, we will give the general form for Wigner functions on the lattice with 
arbitrary lattice size N. As in the previous sections, we assume that the Wigner function 
should satisfy the four conditions (A)~(D) given in Sec. 1. 

It turns out that the Fourier transforms of the Fano operators are much convenient. 

A(n, m) = — J2 e^ (na+m6) A(a, b), (n, m = 0, 1, . . . , N - 1), 

™ a,b=0 
1 N-1 

A(a,b) = j- Yl e-^ (na+m6) A(n,m), (27) 

M n,m=0 

where the site suffixes are suppressed. In terms of these Fourier transforms A, the four 
conditions take simple forms: 

A (n 0) — p^ nai tf 

z - i aia2 \'h u y — c u ai,a2i 
A aiQ2 (0, 77l) 5— ai J,- a 2^rm 

A + (n, m) = A(N -n,N-m), 
Tr(A(n,m) + A(n',m')) = N5 n , n ,5 m , m <. (28) 

We will see there exist infinitely many solutions to the above non-linear equations 
Eqs.([28|). We first show that general solutions can be written in terms of any special solution 
and an arbitrary (N — l) 2 x (N — l) 2 orthogonal matrix. Construction for special solutions 
will be given later. 

^From Eqs. (|28|) we first notice that A(n, 0) and A(0,m) are unique. Denoting a special 
solution by Ao and an arbitrary solution by A, we have 

A aia2 (n,0) = A 0aia2 (n,0) = e l ^ nai 5 ai , a2 , 

A aia2 (0,m) = A 0ai a 2 (0,m) = 5_„ 1+a2 , m . (29) 

These relations and the orthonormal completeness of A imply, for the remaining suf- 
fixes n,m,n' ,m' = 1, 2, . . . , N — 1, each A(n,m) is expressed by a linear combination of 
Ao(n',m')'s. Thus we have 

JV-l 

A(n,m) = Y Unm;n'm'Ao(n',m'), (n,m,n',m' = 1,2, N - 1). (30) 

n',m'=l 

It can be easily shown that A is also a solution to Eqs. (|28|) , if and only if the complex 
(N — l) 2 x (JV — l) 2 matrix U nm - n i m i satisfies the conditions: 

^nm;n'm' = ^N -n,N-m;N-n' ,N -m' j 

U + =W 1 . (31) 

These unfamiliar conditions Eqs.(|3"TD for the matrix U can be expressed in a more familiar 
form by use of an (N— l) 2 x (A^ — l) 2 orthogonal matrix. We introduce an (N — l) 2 x (A^— l) 2 
matrix T as 
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%im;n'm' — 5 n ,N -n> '$m,N -m' , {n, m,Tl ,171 — 1, 2, . . . , N — 1) 



(32) 



with the following properties: 

T = T* = T+ = T-\ (33) 
The conditions Eqs.(|31|) take compact forms as 

U* = TUT, U + = U-\ (34) 
And using the unitary matrix (1 + zT)/y / 2, we write 

where 1Z is an (N — l) 2 x (iV — l) 2 matrix. Now it is not difficult to show that the matrix 
U satisfies the conditions Eqs.(|34]), if and only if the matrix 1Z is real and orthogonal. 

In conclusion, as for the Fourier transforms of the Fano operators with the four properties 
(A)~(D), A(n, 0) and A(0,m) are uniquely given as in Eqs. (p9|) and the remaining other 
ones are given in terms of a special solution Ao and an (N — l) 2 x (N — l) 2 real orthogonal 
matrix 1Z as follows: 

N-l 

A(n,m)= W nm;n / m /A (n',m / ), (n, m, n', m' = 1, 2, . . . , N - 1), 

n',m'=l 

where 

T nm - n ' m ' = 5 n ,N-n'S m ,N-m', (n, TO, Tl', m' = 1, 2, . . . , N - 1). (37) 

Now we will construct special solutions to Eqs.(^8|) and discuss how they are related to 
solutions given before: the solution for odd iV given by Cohendet et. al. and the solution 
for even iV obtained by the reduction method in this paper. 

As we have shown in Eqs.(|28D, A(n, 0) and A(0, m) are unique. They are most conve- 
niently expressed as 

A(n,0) = P n , 

A{0,m) = S m , (38) 

where P is the phase operator and S is the shift operator, which were introduced in the 
previous section: 

P = e i2 ^ ai 8 

Sa\a2 ^01,02 !• (39) 

We can show that the set of operators P n S m (n, m = 0, 1, . . . , N — 1) forms a orthonormal 
complete set as 
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Tr(P n S m ) + (P n 'S m ') = N5 n ,J m , m . (40) 
with the property under hermitian conjugation: 

(P n S m )~^~ = ^w™™- p N ~ n S N ~ m (41) 

The relations Eqs. (|39|) , ( f4~0"|) and fl4~l~D imply that a special solution may be given by the 
form: 

A (n,m) = e ld(n ' m) P n S m , (42) 

if the phase 8(n, m) can be properly chosen so that the condition A + (n, m) = A(N — n,N — 
m) holds. The condition for the phase 6 turns out to be 

8(n, m) = —6(N — n, N — m) + —nm, (mod 2tt). (43) 

We will give some examples. The solution for odd N given by Cohendet et. al. is 
reproduced by the choice: 

. f ^nm (m=even) 
9(n,m) = < %, AT w. r v ) A J . (44) 
I ^ (iv — n) [iy — m) (m=odd) v ' 

After the Fourier inverse transformation, this gives the Fano operator for odd N as 

A aia2 (a,6) = e^ b ^- a ^5 2a , ai+a2 , (45) 

and the Wigner function given in Eq.(0). 

The solution by the reduction method for even N in Sec. 2 corresponds to the choice: 

^nm (n = even, m = even) 

, ^(nm + n) (n = even, m = odd ) fA . 

6(n,m) = { ») \ ' , 46 

' jj(nm — m) (n = odd , m = even) v y 

jj(nm + n — m) [n = odd , m = odd ) 
which gives the Wigner function given by Eq. (p6|) . 



IV. DISCUSSIONS 

We have postulated the four conditions (A)~(D) on the Wigner functions on the lattice. 
After presenting a heuristic reduction method to obtain a special solution on an even site 
lattice, we found the general form of the solutions for any size N , which is parameterized 
by an arbitrary orthogonal matrix. 

Now we will discuss how the quantum mechanical expectation value of a physical quantity 
is related to the Wigner functions. Let us introduce the position, site, operator q and the 
momentum operator ft on the lattice as follows: 

QaiCL2 ^1^01,02) 

i N-1 

Pa,a 2 = jz E &e^(— ). (47) 

iV 6=0 
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From the conditions (A) and (B), it is evident that the expectation value of the physical 
quantity like q n or p m are expressed as 

N-1 

Trpq n = a n W{a,b), 

a,b=0 
N-1 

Tr pp m = b m W{a,b). (48) 

a,b=0 

Differences in many solutions arise, when we evaluate the expectation value for a more 
general quantity involving both q and p. For a general physical quantity, however, it is more 
convenient to use the phase and shift operators P and S instead of q and p: 

P = e*%*, 

S = e^, (49) 

whose matrix elements are given in the previous section. In view of the fmiteness of the 
lattice, the use of P and S is more natural than q and p. Now we calculate the average 
of the classical quantity p n S m with the Wigner function being the weight, the distribution 
function, as 



N-1 

(P n S m ) Wlg ner= £ (e^ a ) " (e^ fc ) ™ W(a, b) . (50) 

a,b=0 



From Eq.([|) and Eqs.(f27|), we find that the right hand side of this equation is nothing but 
the expectation value of the Fourier transform of the Fano operator, and thus we have 

(P n S m ) Wigner = TipA(n,m). (51) 

As we have seen in the previous section, there are infinitely many solutions for A's. If we 
take the form as in Eq. 



A(n,m) = e i9M P n S m i (52) 

then we have 

(P n S m ) Wigner = e ie(n ' m) Tr P P n S m . (53) 

With other solutions for A, we would obtain different operators on the right-hand side. 
Detailed general discussions on this matter, including the continuum limit of the lattice 



Wigner functions, will be given elsewhere [12 



Next we clarify the equivalence between the two ways of site labeling mentioned in Sec. 1. 
In the labeling adopted in this paper, each site is specified by an integer a — 0, 1, . . . ,N — 1 
with the amplitude 0(a). In another way of labeling, which is more appropriate for a system 
of spin s (integer or half-integer) , we label each site by a = —s, — s + 1, . . . , s with the 
amplitude ip{ot). The total number of sites is N = 2s + 1. The correspondence satisfying 
ip(a) ~ 0(a) and ip((3) ~ 0(6) is, up to a constant phase factor, given by 
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; 2vr , 



ip{a) = e l N as (j)(a), a = -s + a, 
#(J3) = e^f'fcb), /3 = -s + b, 

(a,b = 0,l,...,N -1; a,/3 = -s,-s + l,...,s), (54) 

provided that the Fourier transforms ijj and 4> are most naturally defined as 

1 N-l 

m = 7ttE e-^0( a ), (6 = 0, 1, ... , N), 
viv a=0 

f(P) = -J=r E e-**ty(a), 09 = -a, -s + 1, ■ ■ ■ , a). (55) 
yzs + l a= _ s 

Thus it is clear that the Fano operator A s (a, (3) in the spin labeling is related to our A(a, b) 
by the unitary transformation as 

A^ >a >,/?) = e-^(" 1+s ) s A ai+s , a2+s (a + S ,/3 + S )e^( a2+S ) s , 

(a,/3, ai, «2 = — s, — s + 1, . . . , s). (56) 

In fact we can explicitly verify that if A(a, b) satisfies the four conditions, A s (a, (3) is also a 
solution to the four conditions in the spin labeling version, and vice versa. 

However, it should be noted that the Fourier transformation given in Eqs. (|55|) , especially 
in the spin labeling, is not unique. An interesting point in the spin labeling is that the 
Fourier transform Eqs. (|55|) and its inverse transformation naturally lead to the anti-periodic 
boundary conditions for a half-integer spin: ip(a + 2s + 1) = —ip(a) and ip(/3 + 2s + 1) = 
—ip{0). This reminds us of some flavor of the spin statistics theorem as mentioned in [0. 
Discussions on this issue, in conjunction with the reduction method introduced in Sec. 2, 



will be presented elsewhere [[13 
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